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Abstract. We analyze various phases of inflation based on the anomaly- induced effective 
action of gravity (modified Starobinsky model), taking the cosmological constant A 
and k = 0, ±1 topologies into account. The total number of the inflationary e-folds 
may be enormous, but at the last 65 of them the inflation greatly slows down due to 
the contributions of the massive particles. For the supersymmetric particle content, 
the stability of inflation holds from the initial point at the sub-Planck scale until the 
supersymmetry breaks down. After that the universe enters into the unstable regime 
with the eventual transition into the stable FRW-like evolution with small positive 
cosmological constant. It is remarkable, that all this follows automatically, without 
fine-tuning of any sort, independent on the values of A and k. Finally, we consider 
the stability under the metric perturbations during the last 65 e-folds of inflation and 
find that the amplitude of the ones with the wavenumber below a certain cutoff has an 
acceptable range. 
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1 Introduction. 



The importance of inflation becomes more and more apparent in the last years (see, e.g. Q 
for the review). One of the reasons is that the existing and future astrophysical data may give 
a chance to use the early universe as a laboratory for the high energy physics, providing the 
Theoretical Physics a new important information. However, in order to use this opportunity, one 
needs a model of inflation which should be free of any serious phenomenological input. Otherwise, 
the experimental data will mainly deal with this input and not with the underlying high energy 
physics. 

In this paper, we are going to discuss the example of such a model - the new version g, ||, |I| 
of the anomaly-induced inflation (Starobinsky model) ||, ^, |?], [8|, pi This model is based on 
the well known methods of quantum field theory in curved space-time (see, e.g. [11, 12]) and on 
the effective quantum field theory approach. The advantage of the anomaly-induced inflation is 
that it does not contain usual elements of the conventional cosmological phenomenology such as 
inflaton, the precise form of its potential and/or the necessity of a fine-tuning for the initial data. 
The anomaly-induced inflation provides the explanation for both the start of inflation and the 
graceful exit to the post-inflationary FRW-like phase ||. Under a natural supposition about the 
supersymmetry breaking scale we arrive at the magnitude of the Hubble constant which does not 
generate unacceptable amplitudes of the gravitational waves 

The problems of stability play a great role in the anomaly-induced inflation. In particular, the 
transition between the stable and unstable regimes, which could take place due to the breaking 
of supersymmetry || and the quantum vacuum effects of massive fields H, could be responsible 
for the graceful exit to the FRW evolution. In the next sections, we investigate the stability with 
respect to the small perturbations of the conformal factor for various regimes associated to the 
anomaly-induced inflation. One can distinguish three stages of the expansion: 

i) For the almost exponential inflation that takes place at the very beginning (when the effect 
of particle masses is negligible) we consider the spaces with zero, positive and negative space 
curvatures, and with zero, positive or negative cosmological constant. The main part of the results 
can be obtained analytically such that some very important general features are established. 

ii) At the post-inflationary stage, after the supersymmetry breaking, the universe is approaching 
the stable fixed point with a small Hubble constant. In this case the FRW solution is a very good 
approximation for the theory with vacuum quantum corrections, such that these corrections do not 
play an essential role at the later time. We shall verify that the FRW-like solution is an attractor 
for the theory with quantum corrections if the cosmological constant has positive sign. 

Hi) The most complicated is the intermediate transitional epoch between the inflationary and 
the post-inflationary evolutions. This period is characterized by the vacuum quantum effects of 
both massless and massive fields. Indeed, the exact analytic methods for the stability analysis || 
do not apply in this case, and one has to use numerical and approximate analytical methods. 

The paper is organized as follows. In the next section, we present basic concepts of the anomaly- 
induced inflation. In part, we go beyond the results of the previous publications g |, f§. In 
particular, we take into account the effect of the renormalization of the cosmological constant 
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(here, our approach is somehow similar to the one of fl3|| ). In section 3 the stability under the 
small perturbations of the conformal factor is analyzed for the "massless" stage i), mainly through 
the analytical methods. The main new aspect here is that we clarify the role of the cosmological 
constant. In section 4 we verify, using both approximate analytical and numerical methods, that 
the stability of the model with the super symmetric particle content really holds until the scale of 
the supersymmetry breaking. In section 5 we discuss, using the approximate analytical methods 
tested in section 4 and also the numerical integrations, the stability under the metric perturbations. 
Finally, in section 6 we draw our conclusions and also present a general discussion about the possible 
theoretical significance of the model. 



2 Brief survey of the anomaly-induced inflation 

Consider the vacuum quantum effects in the Early Universe, below the Planck scale. Then, the 
appropriate framework is not a string theory but some effective quantum field theory on a classical 
curved background. Let us make some assumptions about a quantum field theory which is valid at 
the sub-Planck energies. It proves useful to concentrate on the asymptotically free theories, that 
helps to avoid the discussion of the higher-loop effects. Thus, in what follows we shall deal with 
the one-loop approximation and only mention possible effects of the higher loops. 

Consider the theory which has Nq real scalars, Nij 2 Dirac spinors, and N\ vectors. One has 
to notice that the vacuum quantum effects originate from the virtual particles. Therefore, the 
numbers Nq, N\/2, N\ correspond to the particle content of the quantum theory, but they do not 
describe the real matter which might fill the Universe. As we shall see later on, the real matter 
has little importance for the anomaly-induced inflation. For simplicity, one can assume that the 
initial universe is empty and that all matter content of the Universe is created during the reheating 
period 

The renormalizable quantum theory of matter fields on curved background includes the action 
of vacuum, and its minimal necessary form is 

Svacuum = SjiD + &EH • (1) 

Here the first term contains higher derivatives of the metric 

S H d = J ' (Pxy/=g {aiC 2 + a 2 E + a 3 V 2 R + a 4 R 2 } (2) 
and the second is the Einstein-Hilbert term 

1 r,4 



Seh = - J d^x^-g (R + 2A) , (3) 

where ai, ..,04, G and A are the parameters of the vacuum action. The renormalization of these 

parameters and the quantum corrections to the action (|l]) are the basics of the inflationary model 

we are discussing here. The discussion in this section will perform in two stages: first we review 

the massless limit and then present a useful approximation Q for the massive case. 

4 Our notations are = diag(-\ ) and 7? M „ = d\Tp V — ... . The sign of the cosmological term in (||) 

corresponds to the positive vacuum energy for positive A. 



3 



2.1 Massless conformal case 



At high energies the masses of the fields are negligible and massless conformal fields serve as 
a good approximation. The vacuum quantum effects correspond to the Feynman diagrams with 
a loop of the quantum matter field and external lines of the classical gravity. Then, in order 
to evaluate the relevance of the mass of the quantum field, one has to compare it to the energy 
corresponding to the external line. In this paper we shall follow Ref. [jnj and associate the energy 
of a quanta of the gravitational field with the magnitude of the Hubble constant H. Therefore, 
the massless approximation is justified for the very large values of H. As we shall see in the next 
subsection, this corresponds to the initial stage of inflation. 

The actions of the free real scalar, Dirac fermion and gauge vector fields are 

S s = ~ J d 4 x^ {(Vc/?) 2 + mV + , (4) 

where (Vip) 2 = g^d^ipd^yp; 

S1/2 = i J d A x^g- { i> [7^ + im f ) $ } ; (5) 



and 
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S 1 = I d A x^{ ~- A F, u F^ \ . (6) 



4 ^ J 

The massless matter fields possess local conformal invariance in the case where scalars couple 
to gravity in a nonminimal conformal way, with £ = 1/6. In what follows, we consider only this 
value, because it provides the traceless stress tensor for the massless theory at the classical level. It 
is well known that at the one-loop level the condition £ = 1/6 holds under renormalization (see e.g. 
fl2|| ), but at higher loop this is not the case [jl4|] . However, even at higher loops the deviation from 
1/6 can be very small []l5|], such that £ = 1/6 can serve as a good approximation. Furthermore, 
massless fermions are conformal invariant, also the gauge fields can be treated as massless and 
conformal at very high energies. 

The renormalization of the Einstein-Hilbert, cosmological and ^—gB? terms in (|l]) is absent 
in the case of conformal massless fields, and from this point of view their presence is not neces- 
sary. Indeed, we have to introduce the Einstein-Hilbert term because it is proved to fit with all 
experimental and observational data at the cosmic scale. For the same reason we introduce the 
cosmological term which was detected in the recent observations [16]. Regarding the J yJ—gR? 
term, we have to accept the presence of this term since it is important for the renormalizability at 



higher loops. But its coefficient may be safely taken very small [[L5| such that it is negligible 
compared to the quantum corrections of the same form (see Eq. (|34j)). For this reason, in what 
follows we shall simply take 04 = 0. Under this choice the action Shd satisfies the Noether identity 

2 6Shd n /»s 
- ^=g t , v - = 0, (7) 

which is usually interpreted as zero trace for the stress tensor of vacuum = 0. At quantum 
level, this identity is violated by the anomaly such that [17] 

2 <5fW 

T =<TH>= -j=9v»-e = ~{wC 2 + bE + cV 2 R). (8) 

V 9 u9fj,v 
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Here is the quantum correction to the classical vacuum action ([!]), the coefficients w, b and 
c are the (3- functions for the parameters 01,02,03 correspondingly 



w 



1 (N JVi /2 Ni 
(4vr) 2 V120 20 10 



(9) 



(No_ 



(4vr) 2 V360 360 



+ 



31iVi 
180 



(10) 



1 (N JV V2 
(4vr) 2 V 180 30 



10/ 



(11) 



Some observation is in order. The direct application of the dimensional regularization jlT]] leads to 
the ambiguous result for the coefficient of the V 2 i?-term in the anomaly. We accept the expression 



for c which emerges from other regularizations such as point-splitting [18] and ^-regularization 



[19]- 



The general solution for the anomaly-induced effective action is j20]. 
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Tind = S c [g^\ + / d 4 x^ { waC 2 + ba(E - -V 2 R) + 2bah A o - 



_l (c + * b)[ R _ 6( va) a -6(VV)] 2 )}, 



(12) 



where 



and 



9flV 9flU 



■ e 



2rr 



-9R 2 



-g[R-6(Va) 2 - 6(V 2 cr)] 2 



A 4 = V 4 + 2 R^V, t V u - ^RV 2 + i (Vi?) V„ 

is the fourth derivative, conformal operator acting on scalars. S c [<7/ii/] = £c[5v] is the "integration 
constant" for the equation (Q). In general, there is no regular method for deriving S c [g^ v ]. For- 
tunately, for some cosmological problems the form of this functional has no importance. Let us 
consider an isotropic and homogeneous metric 



9tw = 9yv ■ a 2 {i) 

where r\ is conformal time and 

ds 2 = g^dx^dx 



drf 



dr^ 



l-kr 2 



r 2 dfl. 



(13) 



(14) 



Then, the conformal functional S c [g^ u ] is constant, for it does not depend on 0(77) and does not 
contribute to the equations of motion. In this particular case, the solution dl2| ) is an exact one-loop 
quantum correction. However, if one takes a more general background, for example 



9nv = 9iw ■ a 2 ! 7 ?) + h^Av, r ) 



(15) 
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where h^ u is a metric perturbation. Then the functional S c [g^ u ] becomes relevant, because 

S c [g^] = S c [g^ v + a~ 2 (r]) • h^(r], r)] 

and the solution ( |l2|) is only an approximation to the unknown general expression. In this and the 
next two sections we shall consider the metrics of the form ( |i~4| ) and in the section 5 the metric 
©• 

The total action with the quantum corrections 

5*4 — S vacuum -\- Tjnd ; (16) 

leads to the following equation for the metric (13): 



o 3 da d 2 / 46 \ da ( 2b\ a Ml fa a 2 k 2A\ . , 

»+ — + ^-( 5+ t)^- 2 H 1 + -)^-^(« + ^ + ^"tJ =»• (17> 

where we used the standard notation for the Planck mass Ml = 1/G. The point over the variable 
designates the derivative with respect to the physical time t, the last is related to the conformal 
time T] by the relation dt = a(r])dr]. 

The equation (|i~7| ) has remarkable particular solutions for all three cases k = 0, ±1: 

k = a(t) = a ■ exp(Ht) (18) 

k = -l a(t) = — sinh(iTt), (19) 

H 

k = +1 a(t) = ^- cosh(ilt) , (20) 

H 

where the Hubble constant H has the form (we do not consider solutions with negative H) 

H=^=.(l±>f^±A) 1/2 . (21) 
V=32i6 I V 3 M p / 



For a negative A, there is only a particular solution with positive sign in (|2l|). 

The solutions for the case of A = have been found by Starobinsky || (see also |^2[ ). According 
to ( |2l| ) and (^), the two solutions with ± exist only for A > 0. Let us remark that in the high 
energy region the sign of A is positive. Let us advocate this point. In the nowadays universe A is 
also positive but very small, due to the cancelation of the induced and vacuum energy densities (see, 



e.g. |23|, [l3j and also 24] for the pedagogical introduction to the formal aspects of the cosmological 
constant problem). 

The induced cosmological constant appears due to the spontaneous symmetry breaking in the 
Standard Model, as a vacuum energy corresponding to the minima of the Higgs potential 

V d = -\m 2 ^ 2 + f -ct>\ (22) 



and has the classical value 

4 

777, 

Knd =< V d >= -— < (23) 

with the magnitude about 55 orders greater than the observable cosmological constant. Obviously, 
the vacuum counterpart must be positive in order the provide a cancelation at low energies. But, 
at high energies (temperatures) above 300 GeV the symmetry in the potential (|22"|) gets restored 
such that the induced part ( f23|) disappears. After that the value of the cosmological constant will 
be given by the vacuum term with the corresponding quantum corrections |TIJ. Similar increase 
of A may be expected in any other phase transition which could occur at high energies. In par- 
ticular, the contribution to the cosmological constant may be expected from the supersymmetry 
breaking in case this breaking is spontaneous. The last observation is that after inflation, when the 
matter content of the Universe gets closer to the equilibrium state, the high-temperature radiation 
dominates over the cosmological term p5fl . 

In what follows, except if this is indicated explicitly, we assume that the cosmological constant 
A in the high energy inflationary epoch satisfies the relation < A <C Mp. Then the expansion of 
( |2l| ) in the small parameter A/Mp gives the two values 



H = J± and H = J-*$--±, (24) 
V 3 V 16vr6 3 ' K J 

where the first solution is exactly that we meet without the quantum corrections. The error in 
both cases has the order of magnitude A/Mp << H. In the present-day universe this means that 
what we observe is really A, without significant quantum corrections. In the inflationary epoch, if 
we compare the second solution (p4|) and the ones in the A = case, the spectrum of H is modified 



by the relatively small A-term such that the two values of H become closer to each other [13] 



As we shall see in the next section, the inflationary solutions (^), (|H]) and ( p0| ) are stable for 



c > and unstable for c < M (see also H ) . According to ( 11 ) , c < is equivalent to the condition 
for a particle content of the theory 



Ni > i N 1/2 + ^ iV . (25) 



The unique way to provide the non-stability of inflation in the case of c > is to introduce 
some additional term or terms into the classical action of vacuum. For instance, the J y/—gR 2 -term 
with the sufficiently large positive coefficient does this job. But, for the first sight, this term is not 
necessary. Let us notice, that the quantum correction (|l^) is valid not only in the Early Universe 
but also at the later stages of the evolution. However, the content No, N1/2, Ni may be different, 
because for the small energies the loops of the fields with large masses decouple. Let us consider the 
present-day universe as an example. The value of the Hubble constant today is extremely small, 
about H rsj 10~ 42 GeV , that is about 30 orders smaller than the mass of the lightest neutrino. 
Therefore, the unique particle which makes a nontrivial contribution to the anomaly is the photon 
0. Then, the relation ( |25| ) testifies a non-stability and the present-day universe is safe without the 
special J ^—gR 2 -term. We shall continue the discussion of this example in the section 3. 
5 If the neutrino would be massless, there could be an equality in (J25|) . 
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The condition fl25| ) enables one to construct a very attractive inflationary scenario H. The 
universe could start in the stable phase, such that inflation starts independent on the initial data. 
The simplest way to provide stability in ( p5| ) is to assume supersymmetry in the high energy region 
and also suppose that the supersymmetry is absent at a lower energy because the sparticles are 
heavy ||. If, in the course of inflation, the magnitude of the Hubble constant decreases, then at 
some point the loops of the sparticles decouple and the matter content Ag. 1/2,1 gets modified. As 
a result, the sign of the inequality (^) changes to the opposite and the universe falls into the 
non-stable inflation with the eventual transition to the FRW evolution. The important problem 
is that why the Hubble parameter should decrease during the inflation? Indeed, the Starobinsky 



solutions (18), fll9|) and (20) are characterized by the constant H. The answer to this question 
is that the value of H really decreases if we go beyond the massless approximation We shall 
consider the effect of massive fields below. 

Another interesting aspect is the possible role of matter. Here we have two different situations. 
The first is the possibility to have a set of heavy massive fields at the beginning of inflation. These 
fields can be remnants of the string phase transition. For our model it is important that this 
massive content does not destroy the stability which holds at the beginning of the stable inflation. 
Let us suppose that these fields create zero pressure. Then the equation (|l7|) is modified by adding 
the term p m a~ s to the r.h.s.. The investigation of the equation with this term and with A = k = 

a a A a z \ c J a° \ c / cr one \a a A J ca° 

has been performed numerically. The result is such that 6VGI1 for Pffi = Mp the system does not 
loose stability and the stabilization performs in a few Planck times. The plots are almost identical 
to the matter less case (see Fig. 2). The reason of this output is that the new term decreases very 
fast during the inflation and becomes irrelevant. 

Second aspect is the transition to the FRW evolution after some of the massive fields decouple 
and the inflation becomes unstable. In order to understand the situation better, let us just sub- 
stitute the corresponding FRW solution a(t) ~ i 2//3 into the equation (26). It is easy to see that 



the "classical" terms - both Einstein and matter, both behave like 1/t 2 , while the quantum correc- 
tion - with higher derivatives, behave like 1/t 4 . Indeed, H(t) — > 0, as usual. Therefore, with the 
properly taken initial data for the unstable inflation, the FRW solution is a perfect approximation 
for the later epochs. In [J|| we have checked that the "proper" initial data emerge naturally in the 
framework of the sharp cut-off approach to the decoupling. 



2.2 Massive case 

Now, let us assume that some of the scalar and fermion fields are massive. At high energies, the 
leading vacuum quantum effect of the masses is the renormalization of the Newton and cosmological 
terms. In Q the renormalization of the Newton term has been taken into account using some special 
approach which gives direct link to the cosmological applications. Here we shall follow the same 
method, but perform the calculations in a more detailed and general form. In particular, we shall 
take into account the renormalization of the cosmological term. 
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The conformal invariance of the actions (|4|) and (||) is violated by the masses and we can not 
use the conformal anomaly to derive quantum corrections. But, this can be changed if we apply 



the conformal description of the massive theory in the framework of the cosmon model [27] (see 
also [^] for other applications of the cosmon model). Let us replace the massive parameters by 
the powers of a new auxiliary scalar field x according to 

2 m 1 i 171 f 

where M is some dimensional parameter. In the original version of the cosmon model J27[| , the 
introduction of the new scalar fields x 1S called to provide the dilatation invariance with the con- 
sequent breaking of the dilatation symmetry. In our case we need a more general, local conformal 
invariance. Therefore, we postulate that the auxiliary field x transforms as 

X^xe~ a , a = a(x), (28) 

independent of the space-time dimension []. After we apply the procedure (|27|), the massive terms 
in (||) and (||) are replaced by Yukawa and quartic scalar interactions between physical fermion 
and scalar fields and the new auxiliary scalar x- The advantage of this procedure is that the new 
actions in both matter and gravitational sectors become conformal invariant. This invariance is 
violated by a conformal anomaly and one can derive the effective action of the background fields 
QauiX using the anomaly-induced effective action scheme. 

The divergences of the theory in the conformal representation have the form 

r 2 = ~(£riJ / d n xV^{wC 2 + bE + cV 2 R+ ^- 2 [RX 2 + Q{dx) 2 \ + ^X 4 }, (29) 

where the coefficients w,b and c are given by Eqs. (||), ( |To| ) and (pd|) correspondingly, and we 
introduced useful notations for the dimensional quantities 

/ = ^X>/ m /' (30) 

and 

Here the sums are taken over all massive fermion / and scalar s fields with the masses m/ and m s 
correspondingly. Nf and N s are multiplicities of fermions and scalars, for example Ni epton = 1 and 
N q Uar k = 3. In order to derive the conformal anomaly, one needs to consider the classical Noether 
identity for the vacuum part of the effective action 

2 faSnac 1 Suae „ , x 

V^g Sgfiv V~9 ox 



6 Other fields transform in a usual way: g M „ — + g M „ exp(2a) , ip — > ip exp( — ^-j^cr) , ip — > ij) exp( — ^— -a). Here we 
use the dimensional regularization, because it helps to simplify the calculations. 
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It is easy to see that, in the theory under discussion, the conformal invariance does not mean 
that the stress tensor is traceless, but instead that the quantity T is zero. Correspondingly, the 
conformal anomaly means < T > ^ instead of usual < T£ > 7^ [|29|J. However, since both 
X-dependent terms in ( p9| ) have the same conformal properties as the square of the Weyl tensor, 
they can be absorbed into the Weyl term such that the problem of deriving anomaly reduces to 
the standard one. Finally, we arrive at the expression 

< T > = - { wC 2 + bE + cV 2 R + X [R X 2 + 6(%) 2 ] + X 4 } . (33) 

In the standard way, using the parametrization g uu = e 2a g U u and \ = e~°"x> one can derive the 
anomaly-induced effective action of the background fields g uu and x 

Tind = S c [g^,x] + J d A x^{wC 2 a + b{E - ~V 2 R)a + 2baAa+ 

+ J f2 [Rx 2 + Hdx)\ + J d'x^gR 2 . (34) 

The last step is to fix the conformal unitary gauge x = X e ~ a = M, such that the classical 
Einstein-Hilbert and cosmological terms acquire their standard form. The one-loop effective action 
becomes 

T = Svacuum + ^ind[X ~> M ] = S HD + S c [g^, M] + 



-g{wC 2 a + b(E - -V 2 R)a + 2baAa} 
3 



3c + 26 
36 



d A x^f=g~R 2 



- / d A Xy 



-ge 



2(7 [i? + 6(W) 2 ] 



1 



16vrG 



j4 / — = 

a xy—ge 



la 



r A 
8^G 



9 ■ cr 



(35) 



It is worth noticing that the expression (|35| ) does not depend on the magnitude of M, but only on 
the fields and parameters of the theory. This effective action differs from the massless counterpart 
( |l6| ) in two respects. First, there are new terms which emerge due to the renormalization of 
the Einstein-Hilbert and cosmological term. Second, the integration constant S c [g^ u ,M] is not 
irrelevant anymore, since we imposed the conformal unitary gauge. Therefore, the expression (f35|) 
can be regarded only as an approximation to the effective action of vacuum for the massive fields 
and we have to learn the limits of validity for this approximation. 

The higher-derivative part of the Eq. ( |35| ) is identical to that for the massless fields, as it 
has to be (see, e.g. p| , |22||). Furthermore, there is a strong link between the anomaly-induced 
effective action Eq.(|35[) and the quantum corrections coming from the renormalization group. The 
expansion of the homogeneous, isotropic universe means a conformal transformation of the met- 
ric 9av{t) = guv exp On the other hand, the renormalization group in curved space-time 
corresponds to the scale transformation of the metric g^ v — > g^ v ■ e~ 2r simultaneously with the 
inverse transformation of all dimensional quantities [O]. For any /j, we have — > fi ■ e T . One can 
compare the a -dependence of the anomaly-induced effective action (|35| ) and the r-dependence of 
the renormalization-group improved classical action 



Si 



vac.imp 



s. 



vacuum 



[P(r)} 



(36) 
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where P = {01,2,3,4 > G, A} denote vacuum parameters of the theory and P(t) = P + (3pT. The 
expression ( |36"1) is a leading- log approximation for the solution of the renormalization group equation 
for the effective action [12] 



r[e- 2T 5a/ g,$ i ,P, A t] =T[g a/3 ,$ i (T),P(T), t j], (37) 

where 5>j are matter fields. It is easy to see that (53) becomes completely equivalent to (^) if 
we set a = const. The coefficient / is a factor of the 0- function for the inverse Newton constant 
1/16ttG and the coefficient g is a factor of the /3-function for the cosmological term A/8ttG. Indeed, 
(|37| ) can be considered as a generalization of the renormalization group improved classical action 

Therefore, we have strong reasons to regard Eq. (|3^) as a leading-log approximation for the 
effective action for the massive fields. This approximation picks up the logarithmic quantum 
corrections and is reliable in the high energy region where masses of the fields are much smaller than 
the Hubble constant H. On the other hand, at low energies the massive fields decouple and their 
quantum effects become negligible. In this paper we shall interpolate between these two regimes 
using the MS scheme and the "sharp cut off" approximation at some specific scale M*. This scale 
is defined such that at H = M* a sufficient number of the sparticles decouple and the inequality 
( |25| ) changes its sign. Within our approximation, the functional S c contains the sub-leading non- 
log corrections and can be disregarded. Let us remark that the M* scale can be different from 
the supersymmetry breaking scale Msusy- In particular, we can suppose that for a GUT model 
Msusy is about 10 16 GeV while the upper bound for M* is 10 14 GeV. This bound appears because 
the amplitude of the generated gravitational waves has an admissible range only for H/Mp < 10 -5 
in the last stage (~ 65 e-folds) of inflation. 

The equation of motion for a(t) has the form 

a + laa + 4<r 2 + 4 (3 - -) aa 2 - 4 - a 4 - 2k (l + — J (a 2 + b)eT 2u - 



Ml 
8vrc 



(<j +2b + ke 



[l- fa)- -fa' 



+ 



M 2 p k 

127TC 



l-ga-g/4) = 0, 



where we introduced useful notations 
~ 16vr/ I ^ V//// 2 



M 



8vrg 

MlA 4?r 



1_ \- iV / m ? . 

3^ Ml ' mp ^ ,„ — „ p __ 



-y 

Ait 



N s m A s 
MlA 



1 N f m ) 

7T 



M 2 p A 



It proves useful to rewrite this equation, also, in terms of a(t) [] 

■ 2 ■■ -2 

1 

7i 



a a a a 
— +3 + 



a a 



4b\ a a / 2b\ a 

5 + — ) - — -2fc(l + — )-= 



(38) 



(39) 



Ml 
8vrc 



- + — + — • l-/-lna 

a a z a z ) 



;,2 



2 a 2 



+ 



MlA 

127TC 1 



g-hxa-g/A) = 0. 



7 We correct a small mistake in the similar expression (without the cosmological constant term) in 
correction does not modify the behavior of a(t). 



(40) 
This 
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The solution of the equation (38) can not be performed analytically in the general form, but it 
is obvious that the approximate solution for small lna(t) can be obtained by the replacement 



M 2 P 



Ml 



Ml 



1-/ In a(t) 



A— >A = A^l-<7lna 



(*)} 



(41) 



in the expression for the Hubble parameter ( |2lD corresponding to the solutions fllq), ( |19| ) and 
(|20|). Another possibility is the numerical integration of the equation (|38|). We present the plots 
produced by the numerical integration at Figure 1 for the case A = 0. These plots were obtained 



using Mathematica program [31] and also FORTRAN and standard algorithms of [3C]. It is easy to 
see that the initial parts of these numerical solutions are identical to the ones which were obtained 
in Q. Let us make one more step in the understanding of these plots and the quality of the 
approximation (^). For the sake of simplicity, we shall consider only the value A = 0. In this 
case one can use (Ell) and easily integrate the equation 



H = H ^l-fa 



Ha 



M P 



<r(0) = 0. 



and find its solution in the form 



a(t) = H t 



Hi 



ft' 



(42) 



(43) 



One might expect that this simple formula should serve as a reasonable approximation only at the 
initial stage when a{t) = lna(i) is small. However, if we compare the plots (a) and (b) at the Figure 
1 and the parabola Fig. lc (flg|), the numerical difference is remarkably small - about 0.1 at the 
maximal point with a max = 10 4 (for the toy model with the MSSM field content and / = 10 -5 ). 
The two solutions are almost identical for all t. Remarkably, this is not a specific feature of this 
particular value of /. If we substitute (|43|) into the equation (|38|) , the result is 



-1.2 • 1(T 8 + (25.9 / - 2.59 • 1(T 8 ) t + (3.14 • 10~ 17 - 3.14 • 1(T 8 /) t 2 + 

+(-1.22 • 1(T 26 + 1.22 • 1(T 17 /) t 3 + (1.48 • 1(T 36 - 1.48 • KT 27 /) t 4 

that is indeed negligible for all situations of interest, when / is very small. This unexpected feature 
of the solution is extremely useful, for we can sometimes use simple expression (fi"3|) instead of the 
complicated numerical solution of the equation ( |38| ) . 

In the last considerations we took into account that / and A are very small quantities. This 
leads to the condition for the magnitude of the quantum corrections to the cosmological constant 



(44) 



Furthermore, one can request that the cosmological constant plays smaller role than the higher 
derivative induced contributions: 



2 E m '- 2 E m /|« M £ WW*) 
/ 



and A < 8vr Mp 1 - / In a(t) 



(45) 



3 The numerical analysis shows that this is a good approximation for A <C / Mp. 
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Indeed, these conditions are satisfied at the beginning of inflation when lna(t) is small and the 
higher derivative terms dominate over the cosmological constant which must be, at most, of the 
order Mi USY - ^ * s convenient to assume the soft supersymmetry breaking scheme such that the 
magnitude of the cosmological constant is much smaller than M^/Mp. When lna(t) becomes large 
and the universe is close to the transition from stable to unstable inflation, we can not rely much on 
Eq.(f|l|) and have to use the numerical methods to verify both stability of the tempered inflationary 
solution and the graceful exit to the FRW phase after the supersymmetry breaks down. 



3 The conditions of stability in the massless case 

In this section we shall consider the stability of the relatively simple massless case (|l7]) . Instead 
of the Eq. ([l?]) one can take the third order 00-component of the equation 

iV - 2 ( R ~ 2A ) 9 ^ = <T ^> ■ ( 46 ) 



Without the matter fields the 00-equation is equivalent to the Eq. ( fl7| ) which we used in the 
previous section || |6|, [32|. The 00-equation can be solved analytically || and the phase diagram 
shows the existence of the single attractor corresponding to the inflationary solution in the stable 
case and various attractors, including H = in the non-stable case with k = A = 0. Unfortunately, 
our equations ( p8| ) or ( |40| ) are much more complicated than in the massless case and there it is 
not possible to solve them exactly. However, there is no real need to do that. At the beginning of 
inflation the stabilization of the exponential solution performs very fast Therefore, at the latest 
stages it is sufficient to check the asymptotic stability of the inflationary solution with respect to 
the small perturbations in the sense of Lyapunov (see, e.g. []33| regarding the stability issues in 
higher derivative models). 

In this section we establish some important general properties of the stability under small 
perturbations and will apply this knowledge for the approximate analytical and numerical analysis 
of the massive case in the next section. 



With respect to the criterion (25) of stability ||, g], there are three relevant questions: i) which 
kind of perturbations one has to consider; ii) whether the choice of k = 0, 1, —1 may affect (f25|); 
iii) whether the presence of the cosmological constant may affect this criterion. 

In order to address the issue i), let us make small perturbations in the equation (|T^) and set, 
for the sake of simplicity, A = and k = 0. First we try the perturbations of the form 

M P 

a(t) = ao(t) + x(t) , where ao(t) = exp(ifi) and H = —= 

V — 167r6 

After simple calculations we arrive at the equation for the perturbations x(t) 

x + 3H'x - (3 + ^) H 2 x - (j + H 3 x + 6(l + ^) H 4 x = . (47) 

Indeed, the factors of H can be easily removed by the simple rescaling of time, and then the 
standard Routh-Hurwitz conditions (see the Appendix) show that the stability can not be achieved 
for any values of b/c. Of course, this contradicts our previous statement (^), and we can conclude 
that the perturbations of a(t) do not serve for our purposes. 
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Fortunately, from the theory of ordinary differential equations (see, e.g. [p4[ ) we know that the 
stability may depend on the choice of the variables. Therefore, let us consider the same problem 
in terms of cr(t). The corresponding equation can be obtained from ( |38| ) by setting A = k = f = 
g = and making small perturbation of the form 

a(t) = ao(t) + y(t) , where cro(t) = Ht . 



The equation for y(t) has the form 



"y + 7Hy 



1(3-* 

c 



H 2 



Ml 
8vrc 



y 



166 TT , MlH 
H H £— 

C 2-7TC 



v = o. 



or, taking the solution for H into account, 



y + IHy + 2(6-- H A y H 6 y 

c/ c 



0. 



(48) 



(49) 



It is easy to see that in this case the stability is achieved when c is positive that immediately gives 
( |25| ) 0. The illusory contradiction between the two results can be clarified if we identify 



y(t) = In 



1 + 



a (i) 



Since ao(i) grows up very fast, the growing modes in x(t) may be negligible because they are 
suppressed by the factor of l/ao(t). Of course, this is the case when c > 0, while for c < there is a 
mode in x(t) which is growing up faster than ao(t). Finally, we learned a very important lesson for 
the consequent study of the complicated massive case: one has to perturb o~(t), but not a(t). The 
calculations for the k ^ and A / cases are a bit more complicated, and we will not bother the 
reader with the details. The condition of stability is always the same as in the simplest situation 
described above. But, it is important that this condition must be the same for k = 0,1,-1 and 
must be independent on A. Let us formulate two simple theorems: 

Theorem 1. The value of A does not influence the asymptotic stability. 

Proof. The equation for the perturbation x{t) or y(t) does not change if the original equation 
for a{t) or o~(t) is multiplied by the factor of a n or exp(ncr). Therefore, we can always isolate the 
cosmological term such that its contribution to the perturbation equation is zero. The dependence 
on A can be also seen in H, but the magnitude of H is irrelevant for stability, because we can 
always renormalize time such that in the new equation H — ► 1. 

Observation. This theorem is valid only for A <C Mp. 

Theorem 2. The choice of k = 0,1,-1 does not influence the asymptotic stability and the 
condition ( p5| ) is valid for all three cases. 

Proof. It is sufficient to apply the corresponding theorem about the asymptotic stability in the 
limit t — > oo (see Appendix). In the equation for the perturbations one can always replace the 
coefficient functions by their asymptotic values. Since both sinh(Ht) and cosh(Ht) behave like an 
exponential at the late times, the conditions of stability for all three cases are the same. 



Of course, the conditions of stability for H(t) are exactly those for a(t). 
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Illustration. The equation for the perturbation in the general massless case is 



'.'/' [t] + 7a y [t] + 



1 26 \ ( .2 .. \ M 2 P 



+ 



+ 



8ff +4 3 



CTn -2k 1 + 



26 



-2(70 



8vrc 



7 ct +8 ( 3 - - ] 



166 . 3 „, A 26 

o-q -4fc 1 + — 

c V c 



4/te~ 2,T0 /i[t] 

2^0 i _ 



Ml 
2nc 



y[t] = o. 



If we neglect those terms in the coefficients of the perturbations which tend to zero at t — > oo, the 
k - dependence disappears. 

Thus, we have learned two lessons from the analysis of the massless case. It is sufficient to 
consider the case A = and k = 0, for this does not modify the stability of the solution. These 
statements can be used in the more realistic theory of inflation based on the effective action for 
the massive fields, but with a proper caution. In fact, the negligible role of A depends on that it is 
a constant. As we shall see in the next section, this is not true for the effective action of massive 
fields. 

Until now, we have considered only the stable inflationary solution with c > 0. In other words, 
we were interested in the high energy region where the supersymmetry is unbroken. Obviously, 
another end of the energy scale also represents a great interest, for it enables one to perform a 
simple and efficient test of the model. Let us consider, again, a nowaday universe with H = Hq ~ 
10~ 42 GeV . According to the recent data [16], this magnitude of the Hubble constant is mainly due 



to the contribution of the cosmological constant. As we have already mentioned in the last section, 
the value of c in the present-day universe is negative due to the contribution of the single massless 
particle - photon. Therefore, we have to check whether the higher derivative terms do not destroy 
the stability of the first solution in Eq. (p4|). Using (fl8|), we arrive at the following characteristic 
equation for the perturbation of H —* H + const ■ exp(At). 



X d + 7HX 2 + 



3- b - 

c 



H 2 



Ml 
8vrc 



A 



c 2nc 



0. 



(50) 



Since the explicit solution of this equation is tedious, let us start from the A 
and the roots of the equation ( |50| ) are 

(o) n xfo) , M P 



A 



0, 



A 



(0) 
2/3 



case. Then H = 



(51) 



At this level, we do not have definite answer to our question about stability. Now, let us look for 



the solution of Eq. (|50| ) making perturbations in (51). Due to the huge difference between H and 
Mp, this approximation is perfect. The result is 



A, 



-AH. 



A 



2/3 



H ± 



M P 



where 



H 



> 0. 



(52) 



It is very nice to see that the cosmological constant A > really stabilizes the solution in the 
low-energy region, exactly as we should optimistically expect. The stability can be verified also for 
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the earlier post-inflationary epoch, when the energy density of vacuum played smaller role than the 
density of radiation (23j. For the later epoch, as we have already noticed, when we substitute 
the FRW solution into Eq. ([17]), the "quantum" terms decrease as 1/a 2 compared to the Einstein- 
Hilbert and matter terms. Therefore, for the later epoch of the expanding post-inflationary universe 
the FRW is a very good approximation. 



4 The stability in the interpolation regime 



At the very beginning of the inflationary period, the evolution goes exactly as in the massless 
case and the condition of stability must be the same (|25j). Later on, when a(t) becomes large 
enough, the deviation from the exponential inflation becomes greater. Finally, the inflation slows 
down such that the Hubble parameter achieves the scale M*, and the massive sparticles decouple. 
Then the condition ( |25| ) breaks down. Our purpose is to see whether the Eq. pH ) remains the 
condition of stability for the modified equations ( |38|) until the point H ~ M*. 

Before starting the numerical analysis, let us perform some analytic investigation of the stability 
in a useful approximation framework. It is sufficient to investigate the behavior of the perturbations 
in the physically important region corresponding to the last 65 e-folds before H achieves the value 
M*. Looking at the plot at Fig. 1, we see that the dependence a(t) is rather smooth. Hence, we 
can divide this region into small intervals (ti,U + Atj) and consider H (t) = Hi = const inside each 
of the intervals. Indeed, H will change from one interval to another, and this may be the source 
of new perturbations. Of course, all this consideration is done for the stable case c > 0, while for 
the non-stable case this approach has no sense. For the purpose of analytic study, we accept the 
approximation (f4l|). Then, the direct calculations give the following equation for the perturbations 
a -> <r + y(t): 



b 'y +hy + b 2 y + b 3 y + b 4 y = . 



(53) 



where 



b = 1 



7H. 



4(3- -)H' 



8vrc 
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2vrc 



4vrc V 1 3 y 



(54) 



and 



M P 
32vr6 



1+1 + 



64vr6 A \i/2 



3 M|/ 



(55) 



If we treat the slowly-varying o~(t) as a constant, then the stability depends on the signs of the 
Routh-Hurwitz determinants (see Appendix). Obviously, D\ = b\ is positive. After some algebra 
we arrive at the expression 



Do 



21M p H 
8vr6 
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77: 



1 + 



64vr6 A \i/2 



3 M 2 p 



> 0. 



(56) 



16 



The general expression for D3 is 



Da 



2\M%H 2 
1Q-K 2 bc 



1 + 



64tt6 A 



+ 1-^ 



1 + 



64tt6 A \i/2 



M 2 p 



4vrc L ^ 3 y J v ; 

Here only the last term may be negative. Let us notice that the values of / and g are small, because 
they are proportional to the ratios of the particle masses and the Planck mass (|39|). Hence, we 
can safely conclude that D3 is positive if the cosmological constant is not too large. Furthermore, 
D4 = 64 • -D3 and the stability of the tempered inflation, within our framework, depends on the sign 
of the same expression which we met in (57) 



H 2 ~f-\l~g. (58) 

The positivity of the last expression is a criterion of stability in the last stage of the interpolation 
regime - just before the decoupling of sparticles and consequent transition to the unstable regime. 
The sign of flag ) depends on the relation between the masses of the fermion and boson components 
and on the way of the supersymmetry breaking. The last defines the magnitude of the cosmological 
constant. If the supersymmetry breaking is spontaneous, the vacuum energy density AM| may 
be as large as Mg USY , that could be, in principle, much greater than M^ . Then, in the case of a 
positive value of g in (|39|), the inflation would be destabilized by the cosmological constant earlier 
than at the M* scale. For the GUT-scale supersymmetry breaking this might be bad, because 
then the magnitude of the Hubble parameter in the last 65 e-folds is too big fj. Thus, for the 
sake of simplicity we shall assume that g is negative or that the supersymmetry breaking is not 
spontaneous such that the vacuum energy density does not exceed M 4 . In each of these two cases 
the transition from the stable to unstable inflation performs, qualitatively, in the same way as for 
A = and the stability of inflation holds until the supersymmetry breaking at the M* scale. 

In the same framework we can investigate the oscillations of the Hubble parameter during 
inflation. This is potentially important issue, since it is related to the possibility of a reheating. 
Let us look again at the equation for the perturbations d53|). Typically, the characteristic equation 
(A4) has one or two couples of complex conjugate solutions 

A = - a ± i/3, (59) 

where a are positive in the stable case. Also, there may be real solutions. For example, in the case 
of the MSSM, A = and / = 10~ 5 , we meet one couple of complex solutions \u 2 ~ — 1.46±3.47i = 
a±i(3 and two real solutions A3 —3.88 and A4 —4.85 • 10~ 6 . The solution for the perturbations 
is 

y(t) = C e~ at cos(/3t + <p) + C x e~ X[it + C 2 e~ Xit . (60) 



10 For the lower scale supersymmetry breaking this does not pose a problem. 



17 



We will be interested in the effect of j3. Of course, for the constant coefficient bk in (|54|) the 
oscillations in the solution are invisible because they are suppressed by the exponentials. However, 
let us remember that the coefficients bk are not constants, for they depend on H = H(t). We can 
consider H as a constant only within a small interval (ti,t% + Atj). When we pass from one such 
interval to another, the value of H changes. The evolution of H will result in the oscillations (|6C|) 
and they may become significant when H changes faster. Let us evaluate this effect. If we define 
bk = H~ l bf i and A = —a + i/3 , — «i/2 to be the solutions of the equation 

6 A 4 + 61 A 3 + b 2 \ 2 + ... + 6 3 A + b 4 = . (61) 

Then ai/2, o. and (3 are constants and moreover 

q 1/2 = Ha 1/2 = H(l - fa)a 1/2 , (3 = H(3 H(l - fa)f3 , a = Hp 9* H{1 - fa)a . (62) 

The last formula shows that the frequency of the oscillations of a(t) and H{t) is greatly decreasing 
during inflation. But, the same concerns also the damping parameter a. Now, let us look at the 
plot Fig. la. We can divide it into three regions with distinct properties: 1) The initial stage 
when masses of the particles is irrelevant. Then, the H (t) is almost constant and the perturbations 
are not visible. 2) The intermediate stage of the "tempered" inflation with a significant effect 
of the particle masses. In this region H(t) changes faster and one can expect oscillations with a 
relevant amplitude. It is easy to see that the amplitude of the oscillations of H (t) has the order of 
mag nitude V" 2 + P 2 C . 3) The final region when a{t) is approaching a plateau and H(t) — > 0, 
where our approximation related to the "sharp cut-off" decoupling obviously fails. In this region 
the sparticles start to decouple and the universe performs the transition to the non-stable inflation 
and eventually into the FRW-like evolution. 

The numerical analysis of stability has been performed using the numerical simulations method, 
and we found that all curves rapidly converge to the non-perturbed one in both the first and the last 
65 e-folds of inflation. We have chosen the initial data <j(ti), <r(tj), <?(ti), cr(ti) for the equation 
( |38| ) using the generator of random numbers independent for each derivative and then integrated 
this equation numerically. The choice of tj has been done in such a way that the corresponding 
Oi = cr(ti) = <jf — 65, where <7j corresponds to H = M*. In all cases we took the same initial value 
a{ti) = 1. The measure of the deviation from the non-perturbed solution H np (t) can be taken, for 
example, in the form 

jj _ I H — Hnp I / H — H n p I ,( H ~ Hn-P \ 

V ^ np J V J \ H-np J 

The result of the time-consuming calculations is that, independent on the choice of the measure for 

the deviation, the solution <to presented at Fig. 1 is stable. Some plots, illustrating the behavior 

of the perturbed solutions are presented at the Fig. 2. The plot of H (t) for the last 65 e-folds at 

Fig. 3 enables one to see the oscillations ^ which are not visible in the scale of Fig. 1 

11 Similar plot with oscillations has been obtained using Mathematica program in the course of preparation of the 
previous article 14]. This plot has not been included into W due to the size limits. 
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Finally, the numerical solutions and approximate analytic investigation of the stability issues 
in this model completely agree between each other. We can observe, at Figure 2, an illustration of 
the stability of the inflationary solution, and at Figure 3 the oscillations of H. We can conclude, 
that our approximate analytic approach is sufficiently robust and is applicable to other problems, 
like the analysis of the gravitational waves and (hopefully) of the density perturbations. 

5 On the amplitude of the gravitational waves 

The most difficult test for a high derivative cosmological model is the behavior of the cosmo- 
logical perturbations. In this section we shall present the preliminary results concerning metric 
perturbations (gravitational waves) during the last 65 e-folds of the anomaly-induced inflation. 
Similar problem was previously studied in J?], |3q| . In both papers, the analysis was restricted to 
the anomaly produced by the massless quantum fields, and the small difference in the equations for 
the perturbations is (in our terms) due to the different choice of the conformal invariant functional 
Scidfiu] m Eq. (|T2|), In the Ref. [^5|, the derivation of equations for the perturbations has been 
based on the anomaly-induced effective action and, in principle, this derivation can be generalized 
for the effective actions Eq. (|34| ) and (|35|). However, this requires a formulation of the local covari- 
ant version of the effective actions (^) and (|35"1), considerable calculations of perturbations for a 
variable H (t) and consequent numerical analysis. At the same time, there is a possibility to evalu- 
ate the stability under the metric perturbations in a more economic way, and that is what we are 
going to discuss here. The detailed study of the gravitational waves will be performed separately. 

Let us make some helpful observations about what we really need to know about the metric 
perturbations. It is well known, that the only phenomenologically important region is the last 65 e- 
folds of inflation. All the gravitational waves which were generated before that are not visible in the 
present-day universe and hence they have no great importance. We shall follow the same procedure 
which proved reliable in the last section when we investigated the perturbations of cr(t). Thus, we 
divide the time interval corresponding to the last 65 e-folds into the small pieces (U,ti + Aij). 
Within each interval (ij,tj + At$), the time derivatives of the Hubble parameter are negligible 
compared to its own magnitude, and we can safely take constant H. In fact, this means that the 
last 65 e-folds can be approximated by the massless model, but with a greatly reduced magnitude 
of the Hubble constant and with much greater magnitude of a = In a q[ As we shall see below, 
the last point leads to the dramatic difference with the original Starobinsky model [M]. 

For the sake of simplicity, let us present the calculations for the A = case only. Using the 
results of |35|], we arrive at the following equations for the perturbations in the massive case 

b h +bi h +b 2 h +b 3 h +b^h+ 

+ n ie ~ 2CT V 2 h +n 2 e- 2fJ V 2 h +n 3 e- 4a V*h = , (63) 

12 At this instant, we have to admit that the choice of the vacuum for the perturbations can be done in the same 
way for the effective actions generated by both massless and massive fields. It would be quite interesting to check its 
validity using some other approach to the derivation of the effective action for massive fields. 
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where h = h(t,x) and 



6 = a\ + w ■ a(t) , b\ = 6H[ai + w ■ &(t)] + 2wH , 

6 2 = UH 2 [ ai + w ■ a(t)] + H 2 (c - 6/2 + 7w) , 

b 3 = 6H 3 [ ai + wa(t)]+H 3 (3c-3b/2 + 5w), 64 = -12H 4 b , (64) 

n\ = — 2H[a\ + w + w ■ <r(t)] , ri2 = — 2[a\w ■ cr(t)] , 713 = a\ + w ■ a(t) . (65) 

We are mainly interested in the dynamics of the perturbations in the last 65 e-folds of inflation 
when the magnitude of a is at least of the order 10 4 . Then, the first observation regarding Eq. (|65|) 
is that the terms with the space derivatives are suppressed by the factor of e 2cr and can be safely 
neglected such that the equation does not depend on the magnitude n of the wavenumber. This 
does not guarantee that the spectrum is flat, because the dependence on n may appear from the 
initial data. For example, in the case of the perturbations of quantum origin the initial data are 
given by ||| 

1 ■ fn ■■ t?I 2 ■■■ 



Even without the n-dependent terms, the equation (|6^) is still complicated, but its analysis 
can be essentially simplified if we remember that a 3> 1. Taking into account that the /3-functions 
(D, ©, (0) all have the same order of magnitude, one can disregard the cr-independent terms in 
d&j]). If we omit also the small term 64, the coefficients of the Eq. ( |63|) become constants and we 
can easily check whether the growing modes are present using the standard methods. Let us make 
a weaker approximation and regard a as a big constant. Then, the Routh-Hurwitz determinants 
are 

D = b , D x = bx, D 2 = 60H 3 [ ai + w a{t)} 2 , 

which are positive for any parameters of the theory; 

D 3 = 360# 6 [ai + w ■ a(t)] 3 + 432H 6 wb[ ai + w ■ a(t)] 2 , 

which is indeed positive for -cr(t) + a\/w > —66/5. Since 6 is a negative number of the order (at 
most) 10, this condition is satisfied if we do not include too large negative coefficient a\ into the 
classical action of vacuum @ . Finally, 

D A = -!2bH 4 ■ D 3 

is positive if D 3 is positive. As we shall see below, for the values n < n c , where n c is a cut-off 
depending on M*, this means that the anomaly-induced inflation is stable under metric pertur- 
bations. This fact is very important for the whole model. If we admit that the transition from 
stable to unstable inflation occurs at H = M* < 10 -5 Mp, then this scale will define the maximal 
amplitude of the gravitational waves and our model fits with the standard restrictions coming from 
the anisotropy of the CMBR. This conclusion is achieved without fine-tuning of the parameters of 
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the model and looks quite robust. In the consideration above we have neglected the derivative H. 
This can be justified, for example, using Eq. (|4l"|). Remember that, by dimensional reason, H must 
be compared to H 2 . A simple calculation gives 



H 2 ' 2(1- fa)' 

For the SUSY GUT the magnitude of / is at most 10 -5 , while the expression (1 — fa) becomes 
comparable to 1CP 5 only at the very end of the stable inflation when the expression (41) and the 
whole approach based on the "sharp-cut-off" needs to be replaced by a detailed description of the 
decoupling mechanism for massive sparticles in curved space-time. Therefore, our approximation 
H/H 2 ~ is reliable in the region where the whole approach is reliable, and can be considered as 
a qualitative hint for the decoupling region. 

The numerical integration of the equation ( |63| ) have been performed for various magnitudes of 
the wavenumber vector n and the initial conditions (|66|). The plots illustrating the stability of 
h(t) are presented at Fig. 4a, where the plots correspond to n < n c = 0.1. As we see, there is 
a very small difference between these plots such that the spectrum is flat. For the larger values 
of n the plots have greater growth at the initial period (see Fig. 4b) After this initial growth the 
perturbations decrease exponentially according to our analytical consideration. But, this decrease 
is too slow and for the large n the waves leave the horizon with the amplitude much greater 
than the initial one. The effect is essentially stronger for smaller values of the parameter /. The 
important observation is that, while the exponential decrease of the perturbations (after the initial 
growth) does not depend on the choice of the vacuum and the approximation H ~ const, the 
dangerous initial growth may change a lot if we take these aspects into account. Therefore, the 
final conclusion concerning the amplitude of the perturbations with large n can be achieved only 
after a more complete analysis. The positive result for n < n c give a hope to meet an acceptable 
behavior of the metric perturbations. 

The main difference with the previous analysis of the perturbations [Q, is that here we have 
to disregard all the u-independent terms in (|6^) . This changes the structure of the equation and the 
growing modes disappear. Let us remark that the numerical analysis of the metric perturbations 
in the first 65 e-folds of inflation shows that, say, the perturbations of quantum origin grow up 
very fast, such that the amplitude increases about 15 orders of magnitude. Indeed, those huge 
waves can not be observed and their existence does not contradict the CMBR data. However, it 
is very important that the magnitude of the perturbations is indeed decreasing with respect to the 
inflationary background metric. As a result, the universe enters into the last 65 e-folds of inflation 
in a very isotropic initial state. All in all, within our framework the anomaly-induced inflation does 
not meet unsolvable problems with the gravitational waves. 



6 Conclusions and discussions 

We have investigated the stability of the anomaly-induced inflation |], [| with respect to the 
small perturbations. The stability holds from the initial stage when the quantum fields may be 
approximately considered massless, until the scale M*, when the most of the sparticles decouple and 
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the inflation becomes unstable. Thus, this model of inflation does not need a fine tuning neither 
for the initial data nor for a graceful exit to the FRW stage. The equations for the gravitational 
waves are essentially different for the first and the last 65 e-folds of the inflation. In particular, 
during the last 65 e-folds the amplitude of the perturbations with n < n c does not grow up and 
the spectrum of these perturbations is flat. Thus, for the choice of the scale M* < 10 -5 Mp, the 
anomaly-induced inflation may be consistent with the CMBR constraints. 

In the period of inflation previous to the transition to the unstable stage, we met a rapid 
oscillations of the curvature, very similar to those which were discovered earlier by Suen and 



Anderson for a similar (technically) higher-derivative model |37| (see also [36|). Therefore, the 
anomaly-induced inflation provides an excellent possibility to describe the reheating during the 
last stage of the stable inflation and in the transitional period. Finally, let us remark that the 
next important steps in the development of this model should be the quantitative description of 
the decoupling process and the corresponding cosmological epoch, and also derivation and analysis 
of the density perturbations. After these two principal problems will be solved, one can start to 
investigate different supersymmetric models and their cosmological consequences. 

At the present level of understanding the model of anomaly-induced inflation looks vert attrac- 
tive from both phenomenological (no fine tuning) and theoretical point of view. Let us discuss 
the last aspect. In the modern high energy theoretical physics there are several drawbacks. One 
example is the supersymmetry and its low-energy breaking. Even if this breaking is a high energy 
effect such that Msusy S> Mp , it looks unnatural that there is not a single couple of the observed 
particles which could be superpartners of each other. This problem becomes even more explicit 
for the low-energy supersymmetry model like MSSM. But, in the anomaly- induced inflation we see 
that exactly this spectrum of supersymmetry is favored by the scheme suggested in |J. Moreover, 
even the "less natural" soft supersymmetry breaking is favorable, because for the spontaneous su- 
persymmetry breaking there could be problems with the decoupling of heavy sparticles. Another 
example is the sign universality of the vacuum /^-functions (^) and (|10|) and the variability of sign 
for ([[l]). As we saw in sections 4 and 5, the conditions b < 0, w > are absolutely necessary 
for the nice stability properties of the inflationary model, while the variability of c is very helpful 
for the successful realization of the transition between stable and unstable regimes |f. It might 
happen that all this indicates to a much more general theory which can give us a correct and 
unified understanding of both vacuum quantum effects and the supersymmetry breaking. From the 
effective quantum field theory point of view, the anomaly-induced inflation shows that there are 
chances to link these different issues in a new way and finally approach this general theory. 
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Appendix. 

This appendix contains a brief summary of the stability issues and especially of the Routh- 



Hurwitz conditions, which were extensively used in the text of the paper. We refer to the book [34| 
for the details. Suppose we know a particular solution xo(t) of the equation 



» 



/(xt"- 1 ),^"- 2 ),...,^) 



(Al) 



In order to check the asymptotic stability of this solution under 



where x = x(t) and x^ 
small perturbations we consider x(t) = xo(t) + y(t) and expand (Al) in the first order in y(t) 



W n) + &iy (n - 1} + 6 2 y (n " 2) + - + 6 n -ii/ + Ky = o, 



(A2) 



where 60 > for definiteness. Let us consider two possibilities: 

i) All hi are constants. Then the problem is reduced to the algebraic characteristic equation 



b X n + 61A + b 2 \ n ~ l + ... + 6 n _iA + b n = , 



(A3) 



because a general solution of (A2) is a linear combination of the terms exp(Aj • t) where Aj are the 
roots of Eq. (A3) and the coefficients are polynomials in t. Hence, if all these roots have negative 
real part 

Re(Ai)<0, (A4) 

the solution yo(t) is asymptotically stable. A useful criterion of (A4) is called the Routh-Hurwitz 
conditions. The conditions are that the determinants Di must be positive, where D\ = b\ and 

/ h 63 65 
b b 2 64 



.Dfc = det 



61 63 
6 b 2 



\ 



hk-i \ 

b 2 k-2 

b 2 k-3 
b 2 k-A 



k > 1 



(AS) 



J 



where k = 2, 3, n and bj = , j > n. For n = 4 we meet four determinants D\ ,2,3,4- 

ii) In the cases k = 1, — 1 of the subsection 3.1 we met the case of the non-constant coefficients 
hi = bi(t). But, there is a general theorem (see e.g. [j34| 1) that the problem of asymptotic stability 
for such case can be reduced to the one with the constant coefficients 6j, provided that there are 
finite limits 6,(t) — ► 6j when t — ► 00. As we saw, this is exactly what happens for k = 1, —1, where 
bi correspond to the k = case. 
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Figure captions 



Figure 1. The plots representing numerical solutions cr(t') (where t' = i/167r) of the 
equation p8| ) for the toy model with / = and MSSM particle content. The plot 

(a) is the same as the one obtained in B. In the region a < M* the two plots (a) and 



(b) are equivalent. The plot (b) is almost identical to the parabola (43). At the top 
point of the plot (b) the value of a equals 1000000, while in case of parabola a equals 
1000001. 



Figure 2. The graphics illustrating the behavior of the perturbation measure D for 
the first (a) and last (c) 65 e- folds of inflation. The plot (b) slightly differs from (a), 
because of the presence of the matter (dust) with the density Mp at the initial point of 
inflation. One can see that the stabilization takes a very short time in all cases. 



Figure 3. The plot of numerical solution for H(t), in the region where the oscillations 
are explicitly visible. 



Figure 4. The samples of the plots of the numerical solutions h(t) in the last 65 e-folds 
of inflation and initial spectrum (|6^) with the normalization h(0) = 1 and / = 10~ 5 . 
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